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Abstract. The structural evolution and dynamics of silver nanodrops Ag2896 (4.4 nm 
in diameter) during rapid cooling conditions has been studied by means of molecular 
dynamics simulations and electronic density of state calculations. The interaction 
of silver atoms is modeled by a tight-binding semiempirical interatomic potential 
proposed by Cleri and Rosato. The pair correlation functions and the pair analysis 
technique is applied to reveal the structural transition in the process of solidification. 
It is shown that Ag nanoparticles evolve into different nanostructures under different 
cooling processes. At a cooling rate of 1.5625 x 10^^Ks~^ the nanoparticles preserve 
an amorphous like structure containing a large amount of 1551 and 1541 pairs which 
correspond to the icosahedral symmetry. For a lower cooling rate (1.5625 x 10^^Ks~^ ), 
the nanoparticles transform into a crystal-like structure consisting mainly of 1421 and 
1422 pairs which correspond to the fee and hep structures, respectively. The variations 
of the electronic density of states for the differently cooled nanoparticles are small but 
in correspondence with the structural changes. 
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1. Introduction 

The research in the field of nanoparticles is the basis for the development of 
nanotechnology [H |2]. This is mainly due to the possibility of modifying the physical 
properties of these systems through the control of the system size. Silver nanoparticles 
are particularly interesting because they have a number of exciting potential applications 
in various fields including electronics and biology [3j. It is also known that the 
structure of a material determines its properties. For instance, in bulk face-centered 
cubic materials the formation of other structures is suppressed kinetically, whereas 
nanoparticles of the same materials exhibit different structural modifications such as 
icosahedron, decahedron and amorphous with a great variety of physical and chemical 
properties [21 HI [5]. Thus, in order to understand the structure of metal nanoparticles, 
obtained from liquid phase, it is important to investigate their structural evolution 
during solidification under different conditions. Molecular dynamics (MD) simulations 
have proved to be among the most effective methods in the investigation (at atomic 
level) of the properties of nanoparticles, which is difficult to carry out experimentally. 
In fact, the structural evolution during cooling of metallic nanoparticles was extensively 
studied employing MD [HI [3, El El [lOl [H] ■ For the case of silver nanoparticles one can 
find some works about, for instance, investigations of the most stable structures and the 
melting process of small silver clusters [12], as well as the superheating of Ag nanowires 
[5]. The freezing of silver clusters and nanowires has also been studied recently |13j . 
In such work the authors conclude that the final structure of Ag clusters of ~2.3 nm 
is a fee polyhedron despite the different cooling rates employed. In contrast, the study 
of the cooling rate dependence of solidification microstructures of silver by Tian et 
al [Hj indicates that the cooling rate plays a crucial effect on the silver structure of 
the solid state. For this case, they find a critical cooling rate for crystal forming of 
r^l.OxlO^^Ks'^ . Baletto et al [15] analyzed the equilibrium structure and melting of 
some magic number Ag nanoparticles, concluding that for relatively large clusters the fee 
polyhedron is the most stable. However, according to experimental results of Reinhard 
et al [16], both icosahedral and fee structures are observed in large (up to 10 nm in 
diameter) free Ag clusters produced in an inert-gas-aggregation source. Hence, details 
of thermal stability, melting and freezing temperature for nanosized non equilibrium 
systems still remain unclear. 

In this work, we investigate the structural transitions of silver nanoparticles during 
fast cooling. The structural changes have been simulated employing a tight-binding 
many-body potential. Furthermore, the possible changes of the electronic properties of 
silver nanoparticles at different temperatures were analyzed by calculating the electronic 
density of states (DOS). The paper is organized as follows. Sec. 2 is devoted to present 
the details of the MD simulations, the structural analysis, and the Hamiltonian model 
for the calculation of the DOS. The results and discussions are presented in Sec. 3. 
Some concluding remarks are provided in section 4. 
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2. Model and methods 

2.1. Potential energy function 

The Ag nanoparticles are simulated by MD methods employing the many-body 
potentials developed by Rosato [IT] on the basis of the second-moment approximation 
to the tight-binding model (SMA-TB). In this framework, the band energy of an atom i 
in a given position is proportional to the square root of the second moment of the local 
density of states. The energy of this atom is then written as a sum of two terms: 

E,,, = j2iE'r''+Er) (1) 

i 

where 

i^^' = -{ E eexp[-2qC-^ - l)]y/\ (2) 

with ^ as an effective hopping integral, rjj- the distance between the atom i and j, r^ 
the cut-off radius for the interaction, tq the first-neighbour distance, and q describes 
the distance dependence of the hopping integral. The second term in eqn. (1) is the 
repulsive energy of the Born-Mayer type: 

K'= E Aexv\-v(-f^-l)l (3) 

The model parameters (^,y4,p, g) are fitted to the bulk properties of the metal. 
The cut-off distance of the atomic interaction is set between the second and third 
neighbour distance. The parameters used in the simulation are taken from [T7]; i.e. 
i = 1.178, A = 0.1028, j9 = 10.928, q = 3.139. 



2.2. Simulation process 

The MD simulations of a silver nanoparticle are carried out for a cubic box without 
periodic boundary conditions, so that the nanoparticle surface is free. The equations of 
motion are integrated in time using a velocity Verlet algorithm. Energy conservation 
with an error less than 1 x 10~^% was achieved with a time step of 6.4 fs. The almost 
spherical nanoparticle was prepared by cutting a spherical region of a desired radius from 
a big face centered cubic crystal. We have considered silver nanoparticles of different 
sizes ranging from 147 atoms up to 2869 atoms, but the results presented here are for 
systems of 2869 atoms, Ag2869 , (in some cases we also present results for nanoparticles 
of 147 atoms, Agi47). It is worth mentioning that 2869 (and 147) belongs to the set of 
magic numbers for icosahedral symmetries [2]. 

In order to obtain an equilibrium liquid-like state for the nanoparticle, we start 
the simulation at 1500 K which is a temperature higher than the equilibrium melting 
temperature of the Ag2869 nanoparticle {Tmeit = 975K). The system is kept at this 
temperature for 10^ time steps (640 ps). Six different quenching processes are then 
carried out. The overall cooling rate was controlled by changing the number of MD 
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steps at each run. As our first cooling process we chose a slow one in which the system 
is cooled from a liquid state at 1500 K to a temperature of 300 K employing 1.2 x 10^ 
MD steps, which corresponds to a cooling rate of 1.5625 x 10^'^Ks~^{kl process). In the 
last one, a fast cooling process, we employ 1.2 x 10^ MD steps, which corresponds to 
a cooling rate of 1.5625 x 10^^Ks~^{k2 process). It is worth mentioning that recently 
Chen et al [10] employed the same cooling rates (/cl, and k2) to study the structure 
and dynamics of a gold nanoparticle of similar size (2112 atoms). Thus, we consider 
these two cooling rates as our extreme cases so that we can compare the results for 
both nanoparticles. The other four cooling rates k were chosen such that kl < k < k2. 
On cooling, the temperature was decreased to room temperature (300 K) by steps of 
AT = lOi^. The internal energy and structural configurations are recorded during the 
simulation. 



2.3. Structural analysis methods 

The pair correlation function. The pair correlation function (PCF) g{r) has been 
widely used to describe the atomic structure in amorphous, liquid and crystalline states. 
This quantity is given by [18], 

^ (n.(r,r + Ar)) 
' pAnr^Ar ' 

where p is the atomic density (N/V), and {ni{r,r + Ar)) is the average number of atoms 
within a spherical shell surrounded by r and r + Ar around an arbitrary atom. 



Pair analysis technique. The common neighbours analysis (CNA) technique, 
introduced by Honeycutt and Anderson [I9] , is a standard tool for the interpretation of 
molecular dynamics simulations of structural transformations. The local environment 
of a pair of atoms is characterized by a set of four indexes {i,j,k,m). The first one 
indicates whether the pair of atoms are closer than a given cutoff distance Tc. In the 
present work Tc is chosen to be equal the semisum of the first and second nearest 
neighbours distance in the perfect fee Ag lattice, that is rc = 0.35 nm. The second 
index j is the number of common neighbours to the two atoms, and the third one is 
the number of bonds between the common neighbours. The fourth index is added to 
provide a unique correspondence between number and diagrams |l9j. For instance, the 
1551 pairs characterize the icosahedral-like local structure whereas the 1421 and 1422 
pairs represent the fcc-like and hcp-like local structures, respectively. Additionally, 
icosahedral and fee systems with structural defects are characterized by the 1541 and 
1431 pairs, respectively. Finally, it is worth mentioning that pairs are named as type I 
if z = 1 and type II otherwise. 
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2.4- Electronic structure 

In order to calculate the electronic properties of silver nanoparticles we employ a tight- 
binding Hamiltonian given by 

R R' 

where is the ortho normal (atomic-like) basis set centered at the site R with one 

orbital s in each site. The matrix elements Hn^jii are defined by 

{6o if R = R 
to if \R-R'\ <r, (6) 
if 1^- ^'1 > re 

with Eo and to as the on-site and hopping terms, respectively, rc is the same parameter 
employed in the CNA technique (see above). The atomic positions {R}, known from 
the MD simulation (see section [273]) . are fixed for these calculations. 

A physically relevant quantity for the study of the electronic properties of the 
system under study is the density of states projected on the atomic position R, also 
known as the local density of states (LDOS). This quantity can be determined from 

nR{e) = lim\--I'mGii^R{6-i'y)\ (7) 

7^0 [ TT J 

where G_r,r(£ — i'j) is the diagonal element of the one-particle Green function G = 
[e — i'-f — H)^^. The recursion procedure is a real-space technique to tridiagonalize 
a symmetric (Hermitian) Hamiltonian such that G can be expressed as a continued 
fraction formed by the so-called recursion coefficients. Details about the recursion 
method can be found in [20] . 

It is worth mentioning that this procedure allows us to collect LDOS in differente 
ways. Considering that an atom in a fcc-solid have 12 first neighbours, for the 
nanoparticle we sum all the LDOS corresponding to atoms with 12 neighbours within 
Tc (fcc-like DOS, FDOS). Thus, in a well ordered nanoparticle with fcc-like structure 
the FDOS will be the main contribution to the total DOS. In this way we can compare 
the structural and electronic changes during the transformation. 

3. Results and discusion 

3.1. Atomic structure of silver nanoparticles 

The nanoparticle Ag2869 is initially in the liquid state forming a nanodrop of 
approximately 4.4 nm diameter. As a result of quenching at different cooling rates 
we obtain, at 300 K, nanoparticles either in a crystal-like (cooling rate kl =1.5625 x 
10^"^ Ks~^) or metallic glass (cooling rate k2 =1.5625 x lO^'^Ks"^) structure, as illustrated 
qualitatively in Fig. [TJ An imperfect nanocrystal with noticeable faceted morphology is 
obtained at a lower cooling rate [kl), as shown in Fig. [D^a) and (b). At a fast cooling 
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Figure 1. Atomic arrangements of the silver nanoparticle Ag2869 at 300 K cooled at 
rates of: a),b) 1.5625 x lO^'^Ks'^, and c) 1.5625 x lO^^Ks-^. The snapshots a) and b) 
correspond to the same configuration but different orientation. Visualization of atomic 
arrangements has been made employing the AtomEye software |2Tj . 

rate {k2) nucleation of equilibrium phases are suppressed, and atoms exhibit random 
arrangements. 

Fig. [2] shows the variation of the average internal energy, Ef, as a function of 
temperature for various cooling rates. At a cooling rate of kl =1.5625 x 10^^i^s~^the 
internal energy undergoes a sharp variation as the temperature decreases from 
approximately 600 to 500 K. As shown in the inset of Fig. [2], the transformation from 
liquid to crystal occurs in a wide temperature range (650 - 450 K): the peak of the 
heat capacity curve dE/dT, located at 527 K, is not sharp. On the other hand, at a 
cooling rate of k2 =1.5625 x lO^^Ks^^the caloric curve Et has no inflexion in the whole 
temperature range. The continuous change of the curve indicates that the system does 
not suffer structural changes; i.e. an amorphous solid is obtained in the undercooling 
state. The reason for this is that high cooling rate restricts the atomic diffusion. This 
behavior is very similar to that reported by Chen et al [10] for a gold nanoparticle of 2112 
atoms. The critical cooling rate, kc-, is the minimum value which is necessary to avoid 
crystallization. For the nanoparticle Ag2869 we estimate kc from the analysis of the time 
evolution of the energy, E{t), as well as its temperature evolution, E{T) (see Fig. [2]), 
and the corresponding PCFs g{r). Thus, we obtain the value of k^. ~ 7.8 x 
Although this value is close to that obtained for solid silver [ll], kc depends on the 
size of the nanoparticle. It is also noticeable that the silver nano glass, obtained at the 
cooling rate of 1.5625 x lO^^fCs"^, is very unstable. For instance, at room temperature 
after 6.4 ps the nano glass transforms to the crystal-like structure in agreement with 

The changes in the atomic distribution within the nanoparticles, during the cooling 
process, is extracted from the pair correlation function. Fig. [3] shows the PCF at several 
temperatures and for the kl and k2 cooling rates. First, in the temperature interval 
from 1500 to 600 K, we note for both cooling processes that the PCFs are identical, 
and reveal the typical structural features of an amorphous liquid with short-range 
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Figure 2. Caloric curves of the Ag2869 nanoparticle for various cooling rates. The 
curve of the fitted derivative dE/dT versus T for kl =1.5625 x lO^^i^rs^-'^is shown in 
the inset. 
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Figure 3. Pair correlation functions of the Ag2869 nanoparticle cooled at a rate of a) 
1.5625 X 10^^Ks-\ b)1.5625 x lO^^Ks-^. 



Atomic and electronic structure transformations 8 



0.20 1 




I I I I \ I I I I I \ I f) on'-' I ' I ' ^ ' I ' I ' ^ 

400 600 800 1000 1200 1400 400 600 800 1000 1200 1400 



Temperature (K) Temperature (K) 

Figure 4. The relative number of bonded pairs 1421, 1422 and 1431 during freezing of 
the Ag2869 nanoparticlc at cooling rate of: a) 1.5625 x lO^^Ks^^, b) 1.5625 x lO^^Ks^^ 
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Figure 5. The Relative number of bonded pairs in the Ag2869 nanoparticle: a) 
1541, b) 1551. kl and k2 are the cooling rates of 1.5625 x 1.5625 x 

10^^ifs~^respectively. 
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topological ordering. At the cooling rate of 1.5625 x 10^^ii's~^the amorphous structure 
of the Ag2869 nanoparticle is conserved down to room temperature. As temperature 
decreases, it is clearly noticeable the splitting of the second peak in the PCF, which 
is characteristic of amorphous structures [22j. On the other hand, the PCF of the 
Ag2869 nanoparticles cooled at a rate of 1.5625 x 10^^Ks~^, below 500 K shows a typical 
structure corresponding to the fee crystal. 

More detailed information on the atomic structure in the nanoparticle can be 
obtained using CNA. At the beginning liquid state (1500 K), the most abundant pairs 
are 1201 (21%), 1311 (21%), 1101 (12%), and 1422 (5.7%), for type I; and 2101, 2211 
pairs, for type II. Figs. H] and [5] show the variation of the relative numbers of several 
principal bonded pairs versus temperature. Quantities are normalized such that the 
total numbers of pairs with i = 1 is unity. For temperatures above 650 K the number 
of principal pairs decreases slowly, but in the same proportion, after increasing the 
temperature which occurs independently of the cooling proceses {kl or k2). This is an 
indication that the structure of the nanoparticle is not changing appreciably, which is 
in agreement with the PCFs analysis (see Fig. [3]). However, for temperatures below 
~625 K the Ag2869 nanoparticle, cooled at a rate of 1.5625 x 10^^ Ks~^, undergo drastic 
structural changes. The number of 1421 pairs increases from 8% to 50% when the 
temperature decreases from 625K to 450 K, while the number of 1431 pairs, which 
represents defect-fcc structure, practically falls to zero. Furthermore, the number of 
1541 and 1551 bonded pairs increases equally, when temperature decreases to 625 K, for 
both kl and k2 cooling processes. Below this temperature the ratio of 1541 and 1551 
pairs falls practically to zero for a relatively slow cooling rate of 1.5625 x while 
for a faster cooling process (e.g. k2), the ratio of these pairs increases, reaching values 
of 12%, and 9% at room temperature, respectively (see figure [5]). It seems that a low 
cooling rate allows the reorganization of the atomic order. As the present simulation 
reveals, the fee crystalline silver nanoparticle of 4.4 nm size is more stable than an 
icosahedral one, which is consistent with the results of Baletto et al [15]. These results 
are in good agreement with the CNA of a gold nanoparticle of 2112 atoms [1^. In the 
case of the small nanoparticle of 147 atoms the number of 1422 pairs, corresponding 
to hcp-like structures, increases up to 35% (300 K) whereas the number of 1421 and 
1551 pairs reach the values of ~ 10% at room temperatures. In summary, for the small 
(147 atoms) and large (2869 atoms) nanoparticles predominate, after the transition, 
the fcc-like and hcp-like structures, respectively. A more detailed study in function of 
the particle size is necessary in order to identify a possible threshold between different 
structures, as was reported in the case of Ni clusters by Qi a/ [6]. 

3.2. Electronic structure of silver nanoparticles 

To determine the DOS of the silver nanoparticles we consider to = leV (reference 
energy), eo/^o = 0, 7/to = 0.2, Ae/to ^ 0.05 (energy step), and 75 (20) recursion 
coefficients for the Ag2869 (Agi47) cluster. 
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Figure 6. total DOS (blue lines), and FDOS (red lines) of the Ag2869 nanoparticle for 
two cooling rates: kl =1.5625 x W^'^Ks~'^{a,) and b)) and k2 =1.5625 x lO^^Ks-^{c) 
and d)). The corresponding temperatures are indicated in the panels. 



Fig. [6] shows the DOS of the Ag2869 nanoparticle for kl=1.5625 x ^and 
k2=1.5625 X W^Ks~^ for both initial (at 1500 K) and final (at 300 K) atomic 
configurations. Comparing the total DOS and the FDOS (fcc-like DOS) we can observe 
that at high temperatures the contribution of FDOS to the total DOS is not remarkable. 
This changes completely at 300 K. For the kl cooling rate, the total DOS at room 
temperature is slightly more similar to the DOS of the fcc-solid sample (not shown 
here), especially for the main peak at ~ — 3i:o, which is not observed in the case for k2. 
This should indicate, in agreement with the analysis of the atomic order (see previous 
sub-section), that the silver nanoparticle with a slow cooling rate becomes more ordered 
(close to its solid counterpart) after decreasing the temperature. The differences in the 
total DOS, at room temperatures, for both cooling rates are more appreciable in the 
small nanoparticle (Agi47), see Fig. [71 For this case the FDOS is different from its 
fcc-solid counterpart. 

4. Conclusions 

The glass formation and crystallization of a supercooled silver nanodrop, 4.4 nm in 
diameter, has been investigated based on the basis of the MD simulation with TB- 
SMA potential. The final structures are highly affected by the cooling rates. Ag2869 
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Figure 7. total DOS (blue lines), and FDOS (red lines) of the Agi47 nanoparticle for 
two cooling rates: kl =1.5625 x m^'^Ks-^{a.) and b)) and k2 =1.5625 x lO^^Ks-^ic) 
and d)). The corresponding temperatures are indicated in the panels. 



nanoparticles, obtained in a supercooled liquid with cooling rates higher than the critical 
value kc = 7.8 x lO^^Ks^^ are very unstable metallic glasses, and nanoparticles resulting 
from relatively slow cooling rates are close packed crystals (fee and hep structures), as 
indicated by the CNA technique. At the cooling rate of 1.5625 x 10^^Ks~^we find a 
discontinuous structural transition near 527 K. The fully crystallized nanoparticle is a 
faceted polycrystal with a number of step consisting of one atomic layer. The study of 
the evolution of the DOS with the temperature indicates that the electronic structure 
depends strongly on the size of the nanoparticle and the cooling rates. 
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